o 

(N 

> 

O 



c3 



i 

c 
o 
o 



(N 
> 

in 

i> 

o 
o 



% 



Density matrix renormalization group study of optical conductivity in the 
one-dimensional Mott insulator Sr 2 CuOs 
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Applying newly developed dynamical density matrix renormalization group techniques at zero and 
finite temperatures to a Hubbard-Holstein model at half-filling, we examine the optical conductivity 
of a typical one-dimensional Mott insulator Sr2Cu03. We find a set of parameters in the Hubbard- 
Holstein model, which can describe optical conductivity for both Mott-gap excitation in the high- 
energy region and phonon-assisted spin excitation in the low-energy region. We also find that 
electron-phonon interaction gives additional broadening in the temperature dependence of the Mott- 
gap excitation. 
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I. INTRODUCTION 

One-dimensional (ID) spin-1/2 quantum spin materi- 
als have attracted great attention because they provide 
typical quantum spin systems with large quantum fluc- 
tuations. A copper oxide compound Sr2Cu03 consists of 
CuO chains with the 3d 9 configuration of Cu ion and the 
2p e configuration of O ion, and it is known as a typical ID 
spin-1/2 Heiscnbcrg system with only nearest-neighbor 
interaction.— Antiferromagnetic exchange interaction J 
between neighboring localized spins is very large with 
J ^0.26 eV (Refs. 2 and 3) but interchain exchange cou- 
pling is very weak. As a result, three-dimensional long- 
range antiferromagnetic order is absent down to very low 
temperature T ^5 K, being regarded as an ideal ID spin- 
1/2 Heisenberg system. 

Sr2Cu03 is a Mott insulator of charge-transfer (CT) 
type. The CT excitation appears at 1.8 eV in optical 
absorption spectral The CT excitation is known to be 
described by a single-band Hubbard model with large on- 
site and nearest-neighbor Coulomb interactions, U and 
V^& via mapping the Zhang-Rice singlet band onto the 
lower Hubbard band in the Hubbard modeli The CT 
excitation gap is, thus, called Mott gap. The optical 
conductivity above the CT excitation gap has been stud- 
ied by using analytical and numerical methods^ - — The 
Hubbard model, however, cannot explain mid-infrared 
absorption inside the Mott gap, which is originated from 
phonon-assisted spin excitation.—^ In order to under- 
stand both the Mott gap and spin excitations on the same 
footing, we need to treat the Hubbard model with opti- 
cal phonons. The simplest model is a Hubbard-Holstein 
model containing Holstein-type coupling of electron to 
the Einstein phonons. 

In this paper, we examine the optical conductiv- 
ity of the Hubbard-Holstein model by using a dy- 
namical density-matrix renormalization group (DMRG) 
methodic combined with a kernel polynomial expansion. 
By treating phonon degree of freedom as a quantum ob- 
ject, we calculate the optical conductivity of the model 
and reproduce both the Mott-gap excitation and phonon- 



assisted spin excitation observed experimentally. We find 
a parameter set describing Si^CuOs. Furthermore, by us- 
ing a low-temperature dynamical DMRG method^ we 
examine the temperature dependence of the Mott-gap ex- 
citation to clarify the effect of optical phonons on spectral 
shape at finite temperature. We find that the presence 
of phonons induces the enhancement of the width of an 
excitonic peak in the optical conductivity. 

This paper is organized as follows. The Hubbard- 
Holstein model in ID is introduced in Sec. II. In Sec. Ill, 
a dynamical DMRG combined with a kernel polynomial 
expansion is explained. The optical conductivity of the 
Hubbard-Holstein model at both zero and finite temper- 
atures are shown in Sec. IV. Summary is given in Sec. V. 



II. HUBBARD-HOLSTEIN MODEL 

A coupling between an electron and breathing phonons 
in cuprates can be mapped onto Holstcin- and Peierls- 
type electron-phonon (EP) interactions^ Since the Hol- 
stein interaction is stronger than the Peierls one, we con- 
sider a ID Hubbard-Holstein model in the present work. 
In addition, we introduce a nearest-neighbor Coulomb 
repulsion that leads to excitonic effects. The Hubbard- 
Holstein Hamiltonian is defined by 

H = -*y^(c|,<7^+i,o- + H.c.) + Uy^nj^nj^ 

i.er i 

+F^(n i -l)(n i+ i-l) 

i 
i i 

where c\ a (ci. CT ) is the creation (annihilation) operator of 

an electron at site i with spin a, and b\ (bi) is the creation 
(annihilation) operator of a phonon at site i. This model 
includes electron hopping, t, on-site and nearest-neighbor 
Coulomb repulsions, U and V, respectively, phonon fre- 
quency, wo, and EP coupling, g. Since the energy scale 



of the optical conductivity is much larger than the dis- 
persion relation of the phonon, we neglect the dispersion 
of the phonon in order to make our discussion simple. 
Keeping cuprates in mind, we take U larger than the 
band width At. 

Here, we briefly mention how to construct the model 
Hamiltonian (fT]) in the case of the cuprates. The lower 
Hubbard band of Eq. (fTJ) corresponds to the Zhang-Rice 
singlet which is derived from a three-band model for Cu 
2>d x 2_ y 2 and O 2p orbitals. Lattice vibration leads to 
changes in the hopping integral of an electron between 
neighboring Cu and O orbitals. The modulation of the 
hopping integral gives rise to a diagonal EP term, i.e., 
Holstein term, in the effective single-band Hamiltonian, 
which is larger than an off-diagonal (Peierls) teri m 13 i 14 

In Hamiltonian (fTJ, we have introduced the U-term, 
V(rii — l)(n i+ i — 1), instead of standard notation, 
VniUi+i. The latter term gives rise to an energy re- 
duction by —2V , instead of — V, when a doubly occupied 
state (doublon) is located at the edge of a chain accom- 
panied by a neighboring unoccupied state (holon). This 
induces a localization of the doublon-holon pair at the 
edge. On the other hand, for the former term, the en- 
ergy reduction is always given by — V for a doublon-holon 
pair irrespective of the location of the pair. Therefore, 
the V(rii — l)(nj+i — 1) term prevents the localization of 
the doublon-holon pair at the edge of a chain, leading to 
the reduction in a boundary effect. 



III. DYNAMICAL DMRG COMBINED WITH 
KERNEL POLYNOMIAL EXPANSION 

We examine the optical conductivity of the ID 
Hubbard- Holstein model [Eq. (fTJ] at half filling. Dynam- 
ical current-current correlation function reads 



XjH = 



i 



v„-*» 



ttN^Z ^ 



Im (n\ j* 



1 



w-H 



ij 



-J \n) , 



(2) 



where N s is the number of electron sites, Z is the 
partition function, j3 is the inverse temperature, j = 
itJ2ia( c i+ia c i-o' ~ H.c.) is the current operator, \n) is 
an eigenstate with eigenvalue e n , and 7 is a infinitesi- 
mally small energy. The optical conductivity is given by 
Xj(ui)/uj. Equation (|2|) for the Holstein- Hubbard model 
is calculated by using a dynamical DMRG method 1 ^ com- 
bined with a kernel polynomial expansion (shown below) 
at zero temperature and a low-temperature dynamical 
DMRG method^ at finite temperatures. 

The dynamical DMRG method employs a multi- 
target procedure. At zero temperature, the multi- 
target states corresponding to the optical conductiv- 
ity [Eq. @] arc |0), j\0), and the correction vector 
[cu — H — e — i'j] j |0), where |0) represents the ground 
state. 



There are several techniques to calculate the correction 
vector. Usually we use a conjugate gradient method ex- 
tended to non-Hcrmitian matrix. In this case, the delta 
function is broaden by Lorentzian with the half width at 
half maximum (HWHM) of 7. Since 7 is usually taken to 
be small but finite number such as 7 = 0.2t in our calcu- 
lation at finite temperature, the tail of the Lorentzian ex- 
tends to a wide range of energy. This tail sometimes cov- 
ers a small spectral weight. Actually, the Mott-gap exci- 
tations at high energy tend to hide the phonon-assisted 
spin excitations at low energy, since the weight of the 
former is more than thousand times larger than that 
of the latter J£ In order to diminish such a undesirable 
tail, we introduce a Gaussian broadening instead of the 
Lorentzian broadening. 

The introduction of the Gaussian broadening is 
achieved by employing a kernel polynomial expansion.— 
Using the Legendrc polynomial Pi(x), we expand the cor- 
rection vector as 



f 



w — H — e,- 

—00 
dx- 



«7 
1 



-J \n) 



00 „i 



oj — x — e n — 17 
Pi(x) 



dx- 



5{x - H)j \n) 

Pi(H s )j\n), (3) 



l l 



where wi = 2/(21 + 1) is a normalization constant for the 
orthogonality of the Legendre polynomial, u/ = w/ju, 
e' n = WH£ n , and H s = whH, wh being a rescaling pa- 
rameter to confine eigenvalues e n within the interval of 
[—1,1]. Note that the Legendre polynomial is defined on 
the real axis within this interval. 

The integration in terms of x in Eq. (J3J) reads 



dx- 



Pi(x) 



'7 



2Q l (u>'-e' n )+iTrP l (Lu'-e' n ),(4) 



where Qi(u) is the Legendre polynomials of the second 
kind. From Eqs. ([3]) and Q, we may calculate the cor- 
rection vector, provided that the Legendre polynomials 
are obtained by a three-term recursive formula. However, 
it is practically impossible to perform integration of I in 
Eq. ([3]) up to infinity Therefore, we need to introduce 
a truncation number L, whose value is usually several 
hundred. 

The truncation of / gives rise to so-called Gibbs oscilla- 
tions in numerical data, which are known to be unphys- 
ical phenomena. Some improved kernel polynomials are 
suggested to remove the Gibbs oscillationsi 16 i 17 In the 
present study, we employ a regulated polynomial^ given 
by a Gauss distribution function, 



(Pl(H.)) s = 



1 



V2ira 2 



dxe 



'Pi(x), (5) 



where a is the HWHM of the Gaussian given by a — 
2tt/L. The regulated polynomial can be recursively cal- 
culated by a three-term recursive formula without direct 



calculation of the Gaussian integral in Eq. ([5]), 

(P l+1 (H 8 )).J\n) = j±±H.{Pln(H B ))-j\ n ) 



xio- 



and 



(Pl +1 (H s )) a j\n) = (2l + l)(P l (H 8 )) a j\ 

+ (PU{H s )) a j\n) 



(6) 



(7) 



where P{(x) = dPi(x)/dx. We note that this regu- 
lated polynomial method gives a Gaussian broadening 
of spectral weight in optical conductivity with HWHM 
of a = a/dJH instead of the Lorcntzian broadening with 
HWHM of 7 in the standard correction vector method. 



IV. RESULTS 

First, we determine a set of parameters of the 
Hubbard-Holstein model [Eq. (JTJ] that reproduces well 
both phonon-assistcd spin excitation and Mott-gap ex- 
citation simultaneously in the optical absorption of 
Sr2Cu03. The phonon energy loq is taken to be ujq = 
0.11 eV from experimental phonon peak at the bottom 
of phonon-assisted spin excitation^ Since an excitonic 
peak seems to exist r- we take a condition for V to gener- 
ate the excitonic peak in ID Mott insulator: V/t = 2&1£ 
Remaining parameters are t. U, and g. From a com- 
parison of optical conductivity between experiment and 
dynamical DMRG calculation, a relation of U/t ~ 8 was 
suggested. 20 ' 21 For g, a diagrammatic Monte Carlo simu- 
lation reported a good description of angle-resolved pho- 
tocmission spectra in two-dimensional cuprates at half- 
filling when g/t ~ 0.4^ though a proper g value may 
depend on the value of ujq. Turning the ratios of U/t 
and g/t together with the value of t, we find a best pa- 
rameter set that can describe both the Mott-gap and 
phonon-assisted spin excitations in different energy re- 
gions. The best parameter set obtained is t = 0.41 eV, 
U = 3.3 eV, V = 0.82 eV, g = 0.16 eV, and u = 0.11 eV. 
The exchange interaction given by a \/U expansion, 
J = U 2 /{U - V) = 0.273 eV, is close to the experi- 
mentally estimated values, J ~ 0.26 eVi 2 ^ We note that 
the parameters of t, U, and V arc different from those in 
Ref. 17, although the difference is not significant. 

Figure [T] shows optical conductivity at zero temper- 
ature for a 24-site chain under open boundary condi- 
tion. We note that similar spectral behaviors in the opti- 
cal conductivity are obtained for a smaller 20-site chain 
(not shown). The Gaussian broadening a is taken to be 
a = 0.04 eV, whose value is enough to smear out discrete 
weights due to finite-size effect [see the inset of Fig.[IJb)]. 




FIG. 1: Optical conductivity in a Hubbard-Holstein chain 
with 24 sites, (a) Phonon-assisted spin excitation, and (b) 
Mott-gap excitation. The broadening factor a = 0.04 eV. 
The DMRG truncation number m = 600 and m = 800 for 
(a) and (b), respectively. Parameter values are shown in the 
panels. 



The number of the states kept in the DMRG procedure 
(m, the DMRG truncation number) is set to be m = 600 
and 800, which is enough to get good convergency. 

Figure [T] (a) shows phonon-assisted spin excitations 
in the optical conductivity. A phonon peak appears 
at ojq = 0.11 eV as expected. Just above the phonon 
peak, a broad structure emerges as phonon-assistcd spin 
excitation. Without EP coupling, we cannot obtain 
this structure. The energy position of a broad peak, 
uj ~ 0.48 eV, is consistent with an experimental valued 
The peak comes from a Van Hove singularity of spinon 
excitations^ 

Figure [T] (b) shows the optical conductivity in the CT 
energy region. A peak appears at oj = 1.75 eV. This 
energy is in agreement with the experimental data^ In 
addition, we can find a hump structure at w = 1.86 eV. It 
is natural to assign the hump structure to a one-phonon 
excitation on top of the ui = 1.75 eV structure because 
of loq — 0.11 eV. Thus EP interaction contributes to the 
broadening of the main peak in the optical conductiv- 

The intensity of the Mott-gap excitation is several 
hundred times larger than that of the phonon-assistcd 
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FIG. 2: (Color online) Temperature dependence of (a) the op- 
tical conductivity and (b) the density of doublon and phonon 
in a 12-site Hubbard-Holstein chain. The truncation number 
m — 400. A Lorentzian broadening factor 7 = 0.2£ in (a). 



spin excitation, which is consistent with the experimental 
data£ Judging from the agreement of calculated spectral 
weights with experimental ones for both the Mott-gap ex- 
citation and phonon-assisted spin excitation, we conclude 
that the Hubbard-Holstein model with the suggested pa- 
rameter set can describe very well the optical properties 
of Sr 2 Cu0 3 . 

Next we examine the effect of temperature on the op- 
tical conductivity in ID Mott insulators. Figure [2] shows 
the temperature dependence of the optical conductivity 
in the Mott-gap energy region for system size N s = 12. 
Since the Herbert space of the N s = 12 system is very 
large due to phonon degrees of freedom, it is impossible 
to perform fully exact diagonalization even though the 
system size is small. Then we use the low-temperature 
dynamical DMRG procedure. In this calculation, we use 
the conjugate gradient method to obtain the correction 
vector with a broadening factor of 7 = 0.2t in Eq. @. We 
employ a parameter set of the Hubbard-Holstein model 
as U/t = 10, V/t = 3, u /t = 0.5, and g/t = 0.25, instead 
of the best parameter set for S^OuC^: U/t = 8, V/t = 2, 
uio/t = 0.27, and g/t = 0.39. We have chosen this new 
parameter set to reduce computational costs because a 
larger uio/t as well as a smaller g/t reduces the number 
of phonons to be included in our numerical calculations. 
The values of parameters are, however, close to the val- 
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FIG. 3: (Color online) (a) The temperature dependence of the 
optical conductivity in a 12-site extended Hubbard chain, (b) 
Comparison between the temperature dependence of the op- 
tical conductivity in the Hubbard-Holstein model (solid lines) 
and that in the extended Hubbard model (dashed lines). 



ues estimated for Sr 2 Cu03. Thus, it is expected that 
the tendency of the temperature dependence of the op- 
tical conductivity is similar to that obtained by the best 
parameter set for Sr2Cu03 used in Fig. 1. The trunca- 
tion number is m = 400. The solid line in Fig. [2] (a) 
shows the result at zero temperature. Since V/t = 3 is 
larger than V/t = 2 in Fig. Q] (b), the excitonic peak 
at uj/t = 5.95 is more isolated from high-energy spectral 
weight consisting of phonon-related states as well as con- 
tinuous states of the Mott-gap excitation. The broken 
and the dotted lines in Fig. [5] (a) represent the results at 
T = O.lt and 0.2t, respectively. We find that the inten- 
sity of the peak decreases with increasing temperature, 
which is consistent with the results obtained experimen- 
tally. 4 It is noted that the peak position shifts toward 
higher u) with increasing temperature. The shift comes 
from the finite size effect, which was confirmed by inves- 
tigating the system size dependence.— 

To investigate the detail of such a temperature 
effect, we calculate the number of phonon per 

site, iV" 1 X)n e_/3e " ( n \ J2i K^i\ n ) an d double occupa- 
tion of electron per site (the number of doublon), 



NZ 1 E„ e~ 0tn (n\ £ 4 4,^A^M- Fi S urc (b) 
shows their temperature dependence up to 0.2£ ~ J/2. 
We find that the phonon number gradually increases with 
increasing temperature, while the doublon number only 
slightly changes. It is natural to consider that phonons 
excited by temperature give some effects on the spectral 
shape of the optical conductivity. 

The intensity of the spectra in both Figs. 2(a) and 3(a) 
decreases with increasing temperature. The reason of 
such temperature dependence comes from the increase in 
the probability of ferromagnetic alignment of neighboring 
spins, which prevents charge transfer process due to the 
electric field. However, phonons may also have a role in 
the temperature dependence of the optical conductivity. 
In order to make clear the role of the phonons, we calcu- 
late the temperature dependence of the optical conduc- 
tivity for the extended Hubbard model without the EP 
interaction in Eq. ([1]) and show the results in Fig. 0(a). 
Comparing Fig. [3] (a) with Fig.[2ta), one can find that the 
peak at finite temperatures in Fig.^a) seems to be wider 
than that in Fig. G3a)- To see this more clearly, we fit 
the peak position of the T — spectrum in the extended 
Hubbard model to that in the Hubbard-Holstein model. 
Here we introduce scaling parameters for energy shift, d, 
and for hight change, h. We obtain d = 0.2 and h = 0.9. 
The rescaled result is shown in Fig. EJb) as dashed lines. 
We find that spectral weight of the Hubbard-Holstein 
model at around uj/t — d = 6.45 , which is higher by u>o 
than the peak position at uj = 5.95, is larger than that 
of the extended Hubbard model. This is due to the pres- 
ence of phonon structure as discussed in Fig. [ljb). The 
same scaling parameters are applied to the T = O.li and 
T = 0.2i spectra for the extended Hubbard model. We 
clearly find that the peak width of the Hubbard-Holstein 
model at T = 0.2t is wider than that of the extended 
Hubbard model. At T = Q.lt, although the difference 
of the peak width between the Hubbard-Holstein model 
and the extended Hubbard model is smaller than that at 
T = Q.2t, the tendency in the peak width is same as that 
at T = 0.2t. This implies that the EP interaction broad- 
ens the peak structure with increasing temperature as a 
consequence of the enhancement of phonon density with 



T as shown in Fig. [2jb). From this result, we find that 
the EP interaction cannot be ignored when we discuss 
the temperature dependence of the optical conductivity 
in ID Mott insulator S^CuOs. 



V. SUMMARY 

We have investigated the optical conductivity in 
ID Mott insulator S^CuOa. We have employed the 
Hubbard-Holstein chain including the Holstein-typc EP 
interaction and Einstein phonon. Using a newly de- 
veloped dynamical DMRG technique combined with a 
kernel polynomial expansion, we have found that, for a 
proper parameter set, our model reproduces simultane- 
ously both Mott-gap excitation at the high-energy region 
and the phonon-assistcd spin excitation at the low-energy 
region in the optical conductivity. We conclude that the 
Hubbard-Holstein model provides a good description of 
Sr2Cu03, and thus the EP interaction plays important 
roles in the electronic structure of Sr 2 Cu03. Using the 
low-temperature dynamical DMRG technique, we have 
found that the EP interaction broadens a peak structure 
in the optical conductivity with increasing temperature. 
This is accompanied by the increase in phonon number. 
Thus the EP interaction also plays an important role in 
the temperature dependence of the optical conductivity 
in the ID Mott insulator such as Sr 2 Cu03. 
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